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Abstract 

Fusion data for 16 0+ 16 are analyzed by coupled-channels calculations. It is shown that the 
calculated cross sections are sensitive to the couplings to the 2 + and 3~ excitation channels even at 
low energies, where these channels are closed. The sensitivity to the ion-ion potential is investigated 
by applying a conventional Woods-Saxon potential and the M3Y+repulsion potential, consisting 
of the M3Y double-folding potential and a repulsive term that simulates the effect of the nuclear 
incompressibility. The best overall fit to the data is obtained with a M3Y+repulsion potential 
which produces a shallow potential in the entrance channel. The stepwise increase in measured 
fusion cross sections at high energies is also consistent with such a shallow potential. The steps 
are correlated with overcoming the barriers for the angular momenta L = 12, 14, 16, and 18. To 
improve the fit to the low-energy data requires a shallower potential and this causes a even stronger 
hindrance of fusion at low energies. It is therefore difficult, based on the existing fusion data, to 
make an accurate extrapolation to energies that are of interest to astrophysics. 

PACS numbers: 24.10.Eq, 25.60.Pj, 25.70.-z 
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I. INTRODUCTION 



A major challenge in nuclear astrophysics is to measure the cross sections for radiative 
capture and fusion reactions with high precision and down to very low energies. Another 
challenge is to develop models that can reproduce the existing data and be used with con- 
fidence to extrapolate the cross sections to the energies that are of interest to astrophysics. 
Examples of reactions where these challenges exist are the low-energy fusion of 12 C+ 12 C, 
12 C+ 16 0, and 16 0+ 16 0. It was recently suggested [1] that the fusion rates that have been 
used in the past for these reactions should be reduced because of a hindrance phenomenon, 
and the implications for stellar burning and nucleosynthesis were investigated in Ref. 
The fusion hindrance has been observed experimentally at extreme subbarrier energies in 
many medium-heavy systems [3| but it has not yet been observed convincingly in lighter 
systems. It is therefore of interest to study the possible evidence for such a phenomenon in 
light systems. 

The theoretical description of light-ion fusion reactions, for example, of the 16 0+ 16 
fusion data [4|], is often limited to optical model calculations. An exception is the study by 
Reinhardetal. fl in which the ion-ion potential for 16 0+ 16 was derived from an adiabatic 
TDHF calculation, and the fusion cross sections that were calculated have served as guidance 
for the extrapolation to low energies. The purpose of this work is to use the coupled-channels 
method to analyze the 16 0+ 16 fusion data of Ref. J4] in order to investigate the influence 
of the couplings to the 2 + and 3~ excitations of the reacting nuclei and to see whether the 
fusion hindrance phenomenon is likely to exist in such a light system. 

The coupled-channels method has been used in numerous analyses of the fusion data for 
medium-heavy nuclei. It has been very successful in many cases in reproducing the data at 
energies near and slightly below the Coulomb barrier, typically down to 0.1 mb or 0.01 mb. 
Thus it has been possible to generate the large enhancement that is needed to fit the data 
at subbarrier energies by including couplings to surface excitation modes, primarily to the 
low-lying 2 + and 3~ states, and to two-phonon and mutual excitations of these states (see, 
for example, the review article Ref. jfj and the proceedings Refs. 3, sj]-) In some cases it is 
necessary to include couplings to transfer channels (see Refs. (3, ls|.) 

The definition of fusion in the coupled-channels approach is usually based on ingoing 
wave boundary conditions (IWBC). They are imposed at a distance somewhere inside the 
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Coulomb barrier, for example, at the minimum of the pocket in the entrance channel poten- 
tial. However, a short-ranged imaginary potential has also been used to simulate the fusion. 
The two ways of defining the fusion give essentially the same result in most cases. This 
supports the view that the fusion is primarily sensitive to the description in the vicinity 
of and outside the Coulomb barrier, and the empirical proximity type Woods-Saxon po- 
tential 9j, which is based on extensive analyses of elastic scattering data and also on the 
M3Y double-folding potential, has served as a very realistic interaction in coupled-channels 
calculations. 

The above view of the fusion process has been challenged by the discovery of the fusion 
hindrance at extreme subbarrier energies. Since the hindrance sets in at a rather high 
excitation energy of the compound nucleus, is was suggested early on that it had to be an 
entrance channel phenomenon [lOj. It was shown in Ref. 11] that the hindrance could be 
explained by adjusting the ion-ion potential at small distances between the reacting nuclei 
so that it produced a shallow potential in the entrance channel. In contrast, the empirical 
interaction of Ref. |9|] produces a relative deep pocket, and the M3Y double-folding potential 
is unphysical because it produces a pocket that is much deeper than the ground state energy 



that the new low-energy fusion data 



of the compound nucleus. It was proposed in Ref. 
offer the opportunity to investigate the radial dependence of the ion-ion potential at short 
distances. 

We have shown that a shallow potential can be constructed by correcting the M3Y double- 
: biding potential with a repulsive term that simulates the effect of the nuclear compressibility 



12]. Thus we were able to explain successfully the fusion data for 64 Ni+ 64 Ni 13(, 28 Si+ 64 Ni 
14j , and 16 O+ 208 Pb 15|, ranging in cross sections from 20 nb and up to 1 b. It is therefore 



of interest to apply this type of potential in the analysis of the 16 0+ 16 fusion data, in order 
to see how it will affect the extrapolation to extreme subbarrier energies. 

With an increased sensitivity (at extreme subbarrier energies) to the potential at small 
distances between the fusing nuclei, one should also be concerned about the validity of 
the basic assumptions in the coupled- channels approach, namely, that the structure input 
is that of the isolated nuclei. For strongly overlapping nuclei, the modes of excitations 
may be better described in terms of excitations of the compound nucleus. Unfortunately, 
it is very difficult in the coupled-channels approach to develop a model that includes a 
realistic transition from a di-nuclear description to a compound nucleus description. Some 
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justification for the coupled- channels method can be found, for example, in calculations 
that are based on the two-center Hartree-Fock method. Thus, in the case of 16 0+ 16 0, Zint 
and Mosel found that the individual shell structures of the two 16 O nuclei survive up to a 



remarkable degree of overlap 



161 ] . Moreover, the ion- ion potential they determined gives a 



very shallow potential in the entrance channel, in qualitative agreement with the empirical 



findings from the analysis of elastic scattering data 171 ]. and also, as we shall see, from the 
analysis of the fusion data. 

II. COUPLED-CHANNELS DESCRIPTION 

The fusion cross sections for the oxygen isotopes are somewhat exceptional in the sense 
that they are rather structureless j^j], whereas fusion data for 12 C+ 12 C [3] and 12 C+ 16 



[19l | contain very rich structures or resonances. Since the coupled-channels calculations 
presented here produce rather smooth and structureless cross sections at low energies, the 
investigations will be restricted to the fusion of 16 0+ 16 0. 

The basic assumptions and ingredients in the coupled-channels description of fusion re- 



actions are summarized below. The structure input for the 2 + and 3~ states in 16 20l. |21| 
is shown in Table I. The channels for the two excited states in 16 become closed when the 
center-of-mass energy is less than the excitation energy, E x fh 6-7 MeV. The asymptotic 
boundary condition for the radial wave function u n i{r) in an open, inelastic channel n with 
angular momentum L is the out-going Coulomb wave, 

u nL {r) oc L {q n r) = G L (r] n , q n r) + iF L (rj n , q n r), for r -> oo, (1) 

where %q n is the asymptotic relative momentum and t] n = Z\Z2e 2 / (hv ) is the Sommerfeld 
parameter for channel n. The out-going wave is here expressed in terms of the regular 
and irregular Coulomb wave functions Fi and Gl, respectively. For a closed channel the 
condition (CQ) is replaced by 

u nL (r) oc W^ VntL+1/2 (2q n r), (2) 

where W- v> l+i/2( z ) is the Whittaker function. 

The coupled equations are solved in the so-called rotating frame approximation. A de- 
tailed discussion of the coupled-equations in this approximation can be found, for example, 
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in Ref. [22j. The boundary conditions that are used at short distances are the ingoing- wave 
boundary conditions (IWBC), which are imposed at the location of the minimum of the 
pocket in the entrance channel, and the fusion cross section is determined by the ingoing 
flux. This model works quite well at energies that are near and below the Coulomb barrier. 

It is a common problem that one cannot always reproduce the fusion data at extreme 
subbarrier energies and at energies far above the Coulomb barrier by using exactly the same 
model assumptions in the two energy regimes. It ha s been suggested that the problem can 

but calculations that demonstrate this 



be solved by considering the effect of decoherence 
point were not carried out. The solution we have used [15( is to supplement the IWBC at 
energies far above the Coulomb barrier with a weak imaginary potential that acts near the 
minimum of the pocket in the entrance channel. The need for such an imaginary potential at 
high energies may reflect the influence of an increasing number of reaction channels, which 
cannot be considered explicitly in a practical calculation. 



A. Standard Woods-Saxon potential 

The real part of the ion-ion potential is commonly parametrized as a Woods-Saxon po- 
tential, 

^ (r) = l + exp[(r-R FCt )/ay (3) 
and the proximity type potential discussed in Ref. [9| (Eq. (40) of Section III.l) will be 
used below. The results of coupled-channels calculations that are based on this potential 



will be compared to data and to calculations that use the M3Y+repulsion potential [12 ]. 
These two types of potentials are basically identical at large radial distances between the 
reacting nuclei, and they produce essentially the same Coulomb barrier height. They differ at 
short distances, as will be shown in the next section. The Woods-Saxon potential produces 
a relatively deep pocket in the entrance channel potential, whereas the M3Y+repulsion 
potential can be adjusted to produce a shallow pocket and a thicker Coulomb barrier. The 

which has been 
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latter two features help explain the hindrance of fusion 
observed in many heavy-ion systems at extreme subbarrier energies. The issue here is 
whether an analysis of the fusion data for 16 0+ 16 will show a sensitivity to the potential 
at short distances. 
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B. Fusion cross sections 

The fusion cross sections for 16 0+ 16 that were measured by Thomas et al. J4] are 
compared in Fig. [TJ to coupled-channels calculations and to the no-coupling limit, i. e., a 
one-dimensional barrier penetration calculation. Both calculations are based on the standard 
proximity type, Woods-Saxon potential 9]. The uncertainty in the data was taken from the 
published figures, except at center-of-mass energies larger than 8 MeV, where an (arbitrary) 



uncertainty of 5% was adopted, because it was not possible to read the sma 



uncertainty. Earlier measurements do exist, see for example Refs. 
not be shown here. 



24 
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1 experimental 



26| , but they will 



in which the radius has been adjusted to 
, has the parameters: Vq = -42.14 MeV, 



The proximity type, Woods-Saxon potential 
provide the best fit to the 16 0+ 16 fusion data 
R po t= 6.083 fm, and a = 0.602 fm. The solid curve (CCC) in Fig. [TJ is the coupled-channels 
result one obtains with this potential and it has a x 2 /N = 1-5- The dashed curve shows the 
no-coupling limit (NOC) one obtains with the same potential. The fit to the data in the 
no-coupling limit can be improved by adjusting the radius of the potential. The best fit is 
achieved for R pot = 6.133 fm and has a x 2 /N = 4.2. This is much larger than the x 2 /N 
= 1.5 obtained in the coupled-channels calculation and shows that the couplings to the 2 + 
and 3~ states do play a significant role. 

It may be seen in Fig. [TJ that the measured cross sections fall off faster with decreasing 
energy than predicted by the coupled-channels calculation at the lowest energies. This is a 
signature of the onset of the fusion hindrance phenomenon discussed earlier. The hindrance 
will be explored further in the following sections. 



III. THE M3Y+REPULSION POTENTIAL 

The calculation of the M3Y+repulsion potential is described in Ref. [3] but some of 
the essential features are summarized here. First one calculates the M3Y double-folding 
potential (including the exchange term). This requires as input the proton and neutron 
densities of projectile and target. The densities of protons and neutrons were assumed to 
be identical for 16 O and the form of density was assumed to be a fermi function with radius 
R=2.5 fm and diffuseness a=0.52 fm. The radius was adjusted so that the measured RMS 
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charge radius of 2.737(8) fm [27] was reproduced. 

The repulsive term associated with the nuclear incompressibility is also obtained from 
the double-folding procedure using a repulsive effective NN interaction of the form V Tep 8(r). 
The density that is used in connection with the repulsive term has the same radius as the 
ordinary densities mentioned above but the diffuseness a rep is chosen differently. Thus there 
are two parameters in the calculation of the repulsive interaction, the stren-h^ and 
the diffuseness a rep of the density. They are constrained, as explained in Ref. [12], so that 
the total nuclear interaction for completely overlapping nuclei, U^{r = 0), is consistent the 
equation of state e(p) at normal nuclear matter density p, and e(2p) at twice the nuclear 
matter density. This condition was expressed in Ref. [12j by the relation 

U N (r = 0) = 2A a (e(2p) - e{p)) « ^ K. (4) 

Here A a is the mass number of the smaller nucleus, so the equation expresses the change in 
energy one has by embedding the smaller nucleus inside the larger. The last approximation 
relates this change in energy to the nuclear incompressibility, K = 9p 2 (d 2 e(p)/dp 2 ). 

The entrance channel potential obtained from the M3Y+repulsion potential is illustrated 
in Fig. [2] for a range of values of the diffuseness parameter a rep . The strength of the repulsive 
interaction V Tep was adjusted in each case to produce the nuclear incompressibility K = 234 
MeV. This is the value that has been obtained from the Thomas-Fermi equation of state for 
symmetric nuclear matter 28(. The smallest value of a rep , which is 0.3 fm, produces a pocket 



that is as deep as the energy of the compound nucleus 32 S. The largest value, a rep =0.41 fm, 
produces a pocket at 2.4 MeV and has a Coulomb barrier of 10.01 MeV. 

The thin solid curve in Fig. [2] is the entrance channel potential one obtains with the pure 
M3Y potential (including the exchange term). It has an unrealistic and extremely deep 
pocket, which is far below the energy of the compound nucleus. The entrance potential, 
which is based on the Woods-Saxon potential discussed in the previous section, is shown by 
the lower thick dashed curve. It is slightly deeper than the ground state of the compound 
nucleus and it has a Coulomb barrier of 10.10 MeV. 

Finally, the upper thick dashed curve in Fig. [2] is the shallow Gobbi potential 



which, by the way, is in surprisingly good a 
two-center Hartree-Fock calculation of Ref. 



greement with the potential obtained in the 



161 ] . It is seen that the M3Y+repulsion and 



Gobbi potentials have almost the same depth but the thickness of the Coulomb barrier is 
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different. It turns out that the thicker barrier provided by the M3Y+repulsion potential 
gives a much better fit to the fusion data when applied in the coupled-channels calculations. 



A. Fusion cross sections 

The thick solid curve in Fig. [2, which is based on the diffuseness parameter a rep = 0.41 
fm, is the entrance channel potential that provides the best fit to the fusion data in the 
coupled-channels calculations. The calculated cross section is shown by the solid curve in 
Fig. [3] and it has a ^/N = 1.3. The quality of the fit is only slightly better than what 
was obtained in the previous section using the standard Woods-Saxon potential. The slight 
improvement is difficult to see but it is achieved mainly at the lowest energies. 

The y 2 /N is shown in Fig. H] as function of the diffuseness parameter a rep . There are 
two minima, one at a small value, a rep ~ 0.325, and one at a rep = 0.41 fm, which is by far 
the best solution. It is of interest to compare a rep to the values that have been used for 
other systems. Thus for the 64 Ni+ 64 Ni system we obtained the best fit to the data for a rep 

n _ „ D _ . . n . „ , 

:or the very asymmetric 



15| 



= 0.403 fm (12| . For 28 Si+ 64 Ni the value was 0.392 fm [14j, and 
system 16 O+ 208 Pb we had to use the smaller value a rep = 0.35 fm 

The discrepancy with the data is emphasized in Fig. [5] where ratios of the measured 
and calculated fusion cross sections are shown. The solid circles are the coupled-channels 
results and the open circles show the results in the no-coupling limit obtained with the 
same potential, namely, the M3Y+repulsion potential with a rep = 0.41 fm. The effect of 
the couplings to the 2 + and 3~ excitations is to bring the cross section ratio closer to one. 
However, there are still some minor deviations from one. For example, the ratio of the 
measurement and the coupled-channels calculation (solid circles) shows a decreasing trend 
with decreasing energy below 8 MeV and it is less than one at the lowest energy point. This 
is a signature of the experimental fusion hindrance with respect to the coupled-channels 
calculation. 



IV. THE S FACTOR AT LOW ENERGIES 

It is unfortunate that the quality of the fit of the coupled-channels calculations to the 
16 0+ 16 fusion data is essentially the same whether we use the Woods-Saxon or the 
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M3Y+repulsion potentials. From the empirical knowledge of the fusion hindrance phe- 
nomenon [2I one would have expected that the M3Y+repulsion potential would provide 
a much better description of the low-energy data. However, the improvement is modest. 
One would need measurements at even lower energies in order to be able to see a stronger 
sensitivity to the ion-ion potential at short distances. 

A good way to emphasize the low-energy behavior of the measured and calculated fusion 
cross sections is to plot the S factor for fusion defined by 

S = Ecm.tr f exp(2vrr]), (5) 

where 77 is the Sommerfeld parameter. The experimental S factors are compared in Fig. [6] to 
the two coupled-channels calculations that were discussed earlier. The top dashed curve is 
based on the Woods-Saxon potential, whereas the solid curve is based on the M3Y+repulsion 
potential. The latter provides a slightly better fit to the data at the lowest energies but the 
error bars are large so the overall improvement in terms of a x*/N is modest, as discussed 
in the previous section. 

The two coupled-channels calculations shown in Fig. [6] start to deviate as the energy is 
reduced. The calculation that is based on the M3Y+repulsion potential (the solid curve) 
develops a maximum near 4 MeV. The reason is that the entrance channel potential has a 
pocket at 2.4 MeV, and this forces the S factor to vanish below that energy when the fusion 
is determined by IWBC It is interesting that the value of the S factor at 4 MeV (solid curve) 
is in fair agreement with the prediction of the adiabatic TDHF calculation |5j. However, 
the S factor for fusion obtained in the adiabatic TDHF calculation keeps increasing with 

n 

decreasing energy p]. 

It is not clear a priori whether the S factor for the fusion of 16 0+ 16 should develop a 
maximum at low energy. It does not have to do that because the ground state Q value for 
producing 32 S is positive. It is only for negative Q-values one can argue that the S factor 
must have a maximum at some positive center-of-mass energy [3J]. 

There is an alternative extrapolation method [s| which is based on the logarithmic deriva- 
tive of the energy-weighted fusion cross section, 

This quantity has a nearly linear dependence on energy at extreme subbarrier energies 
in most of the medium-heavy systems that have been studied experimentally. The linear 



dependence makes it fairly easy to extrapolate the data to the energy where the S factor 
has a maximum j^]. 

A better parametrization of L(E) was adopted in Ref. lL By considering all of the 
50 data points that have been measured below 8.5 MeV {4, 24, 25, 26] it was concluded 
that the S factor for the fusion of 16 0+ 16 must have a maximum close to 7 MeV. The 
extrapolation to lower energies that was obtained in Ref. jl| is shown in Fig. [6] in terms of 
the S factor by the lowest, thick curve. The low energy fusion cross sections predicted by 
this extrapolation are even more suppressed than the coupled-channels calculation that is 
based on the M3Y+repulsion potential. 

Apparently, there are certain features of the data that are not reproduced by the coupled- 
channels calculations presented here. Some indications of that can be seen in the cross 
section ratios shown in Fig. [5j For example, the seven lowest data points form an isolated 
group which is disconnected from the rest above 8.5 MeV. It is not clear what causes the 
discontinuity; is it a remnant of a resonance or is it an experimental problem? In any 
case, one can adjust the M3Y+repulsion potential so that the coupled-channels calculations 
reproduce the the energy dependence of the seven lowest points, i. e., so that the cross section 
ratio <Jexp/&caic becomes a constant. This can be achieved with a diffuseness parameter in 
the range a rep = 0.425 - 0.43 fm, which produces a pocket in the entrance channel potential 
in the range of 4.5 to 5.2 MeV. The S factors obtained from such calculations are also shown 
in Fig. [6j The results are in fair agreement with the extrapolation method proposed by 
Jiang et al. [1]. That is not surprising because the latter extrapolation was also based on 
low-energy data. 

Fusion should in principle be allowed down to zero energy because the ground state of 
the compound nucleus 32 S is at a much lower energy (see Fig. [2]). To describe the fusion at 
such low energies would require an extension of the model used here, for example along the 
lines proposed in Ref. [291 ] . 



V. FUSION AT HIGH ENERGIES 



Another way to test the ion-ion potential is to compare to cross sections that have been 
measured at energies far above the Coulomb barrier. Here the data are often suppressed 
compared to calculations that are based on a conventional Woods-Saxon potential, with a 
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relatively deep entrance potential [30j. We have previously shown that the shallow entrance 
channel potential, produced by the M3Y+repulsion potential, gives a better description of 
the high energy fusion data for 16 O+ 208 Pb 15) . However, it was necessary to supplement 
the nuclear interaction with a weak, short-ranged imaginary potential. It is interesting that 
the same conclusions apply to to the high-energy fusion of 16 0+ 16 0. 

The results of coupled-channels calculations that are based on the M3Y+repulsion po- 
tential and a short-ranged imaginary potential that acts near the minimum of the potential 
pocket are compared in Fig. [7] to the data of Tserruya et al. 31]. The data exhibit an 
oscillatory or step-wise increasing behavior which was also seen in the earlier data by Kolata 



et al. 



32[| This behavior is qualitatively reproduced by the calculation (solid curve). The 



step-wise increase in the calculated cross section is correlated with overcoming the potential 
barriers for L = 12, 14, 16, and 18. This can be seen by comparing to the thin dashed curves 
which show the cross sections one obtains by imposing different values of the maximum an- 
gular momentum L max in the calculations. Note that only even values of L are considered 
for a symmetric system. 

The coupled-channels calculations that are based on the deep Woods-Saxon potential are 
shown by the upper thick dashed curve in Fig. [7J It is seen that the data are suppressed 
compared to this calculation and that the step- wise behavior sets in at a higher energy and 
a higher angular momentum (L=16 to be precise, compared to L=12 in the solid curve.) 
Thus the high energy 16 0+ 16 fusion data show a clear preference for the shallow potential 
produced by the M3Y+repulsion interaction. 

The analysis of the elastic scattering data for 16 0+ 16 by Gobbi et al. [17i] also revealed 
the need for a shallow potential. The potential they obtained is illustrated by the upper 
thick dashed curve in Fig. [2j The minimum of the pocket is in this case at 0.78 MeV, which 
is slightly deeper that the 2.4 MeV pocket produced by the M3Y+repulsion potential (with 
o re p = 0.41 fm.) Thus it appears that both the elastic scattering data and the high energy 
fusion data prefer a shallow pocket in the entrance channel. 

It should be mentioned that the Gobbi potential does not provide a good description of the 
low-energy fusion data by Thomas et al. [4j, although it has a shallow pocket. The reason is 
that the Coulomb barrier is not as thick as the one produced by the M3Y+repulsion potential 
(see Fig. [2j) As a consequence, the fusion data are hindered compared to calculations that 
are based on the Gobbi potential. 
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VI. CONCLUSIONS 



It has been shown that the calculated fusion cross sections for 16 0+ 16 are sensitive to 
couplings to the 2 + and 3~ excited states of 16 even at low energies, where the excitation 
channels are closed. Unfortunately, the overall quality of the fit to the fusion data by Thomas 
et al. Q| is not very sensitive to the ion-ion potential at short distances between the reacting 
nuclei. It is only at the very lowest energies that there is a preference for a shallow potential 
in the entrance channel. 

The potential that gives the best fit to the fusion data by Thomas et al. is the M3Y 
double-folding potential which has been corrected for the effect of the nuclear incompress- 
ibility. This M3Y potential is calculated with a density that is consistent with the measured 
charge radius of 16 0, and it produces a very realistic height of the Coulomb barrier. The 
repulsive interaction that simulates the effect of the nuclear incompressibility is calculated 
with parameters (the nuclear incompressibility and a diffuseness parameter) that are sim- 
ilar to those that have been used previously to reproduce the low-energy fusion data for 
medium-heavy systems. 

The fusion cross sections obtained in coupled-channels calculations are in fairly good 
agreement with measurements at high energies when the calculations are based on the shal- 
low M3Y+repulsion potential. In particular, the oscillatory or step-wise increasing behavior 
is reproduced very well, whereas the calculations that are based on the deeper Woods-Saxon 
do not reproduce the data. The evidence for a shallow entrance channel potential is corrob- 
orated by the empirical optical potential for the elastic scattering obtained by Gobbi et al. 

121- 

The S factor obtained in the coupled-channels calculations that give the best fit to the 
data by Thomas et al. has a maximum at a center-of-mass energy near 4 MeV. The value of 
the S factor at this energy is close to the value that was predicted more than 20 years ago in 
an adiabatic TDHF calculation. However, if the potential is adjusted to improve the fit only 
to the low-energy data, one obtains an even stronger hindrance of fusion at lower energies. 
This is in qualitative agreement with the empirical extrapolation proposed recently by Jiang 
et al. To confirm the hindrance experimentally one would have to measure the fusion 
cross section down to an energy of 5-6 MeV. 
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FIG. 1: The fusion cross sections for 16 0+ 16 measured by Thomas et al. are compared to 
coupled-channels calculations (CCC) and to the no-coupling limit (NOC). Both calculations are 
based on the Woods-Saxon potential described in the text. 
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FIG. 2: The M3Y +r e P nl S ion —channel potential is S h„ W n for * =0,! f m , and lot ^ _ 
0.3-0.4 fm in steps of 0.025 fm. The Woods-Saxon, the Gobbi [13], and the pure M3Y entrance 
channel potentials are also shown, and the energy of the compound nucleus 32 S is indicated. 
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FIG. 3: The fusion data for ib O+ ib O 
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J] are compared to coupled-channels calculations (CCC) 



and to the no-coupling limit (NOC). Both calculations are based on the M3Y+repulsion potential 



with a rep = 0.41 fm. 
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FIG. 4: The x 2 /N obtained from the 16 0+ 16 fusion data ^ and coupled-channels calculations. 
The x 2 /N is shown as function of the diffuseness parameter a rep , which determines the repulsive 
term in the M3Y+repulsion potential. 
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FIG. 5: Ratio of experimental 4] and calculated cross sections shown in Fig. [3j 
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FIG. 6: S factors for the coupled-channels calculations shown in Figs. [Hand [3] are compared to 
the data [4J. Also shown are coupled-channels calculations that are based on shallower potentials, 



with a 



rep 



0.425 and 0.43 fm, respectively. The diamond is t 



lowest ( hick curve is (lie extrapolation made by Jiang et; al. [l] 



le adiabatic TDHF results [5|. The 
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FIG. 7: The measured fusion cross sections for 16 0+ 16 (solid circles 4|, diamonds: 
coupled-channels calculations that are based on the Woods-Saxon (top dashed) and M3Y+repulsion 
potential (with a rcp = 0.41, solid curve). The thin dashed curves show the dependence on the 



maximum angular momentum L max . All calculations include a short-range imaginary potential. 
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